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9^' Abstract 

CN . In the paper are considered stationary (Bloch) states of a particle, in the 



field of periodic biparabolic type potential. It is shown that while the 
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in contrary to the expected decreasing. This "anomalous" behavior is more 
pronounced for the near-top bands and monotonically decreases for the higher 
or lower ones. 

I. INTRODUCTION 

Bloch functions are the basis of periodic system theory, particularly of crystalline solid 
physics and optical lattices 0. For a high-precision approximation of a standing wave 
(sinusoidal potential), by one of the authors was suggested a biparabolic form for the po- 
tential 0, Bloch functions of which are expressed by confluent hypergeometric functions. 
The approximation for the potential has been used to calculate the spontaneous emission in 
the field of resonant standing wave, the presence of light-induced anisotropy in that emis- 
sion was shown |Q. Also it was used for calculating the temperature of ideal gas BEG in a 
periodic field and was shown that the critical temperature decreases with the deepening of 
the potential 0. 



In this paper the biparabohc periodic potential is used to investigate the Bloch states 
themselves, more concretely, for identifying the behavior of modulo wave function as a 
function of energy in the limits of a specific energy band. Specially we are interested in 
so-called " underbarrier" (classically forbidden) states. 

The way of dependance of underbarrier states on energy is well-known for a single po- 
tential barrier and for transitions from one energy band to the neighboring ones: the under- 
barrier wave function decreases conjointly with the decreasing of energy of a particle. But 
what is the kind of dependance in the limits of a sole (definite) energy band? Surprisingly 
this question has slipped from investigators sight; at least in scientific literature available 
to us we did not find works concerning this question. As a formal argument for examining 
this question, for us appears the fact that Bloch states are generated through the interfer- 
ence of particle de Broglie waves transmitted and deflected from many (infinite in the limit) 
potential barriers, and hence interference, besides the energy, also will play an essential role 
in establishment of regularities. Of course, the deflection and interference of underbarrier 
states take place in case of a single-barrier potential too, but they are small there and as 
a dominating behavior will appear the energetic behavior. Consequently, the inter-band 
regularity is not under necessity to be repeated for intra-band states too. Really, in case of 
periodic potential, when from different barriers multiple-wave interference takes place, as a 
dominating one can appear the interference behavior , and the dominant role of interference 
will be expected for intra-band regularities. Than, if the interference behavior does not 
coincide with energy one, the general behavior will mimic the interference but not energy 
regularity. Our calculations show that indeed, when particle's energy decreases in limits 
of a deflnite band, the modulus of the wave function of underbarrier (classically forbidden) 
region does not decrease, as should be expected from energy reasonings, but increases. The 
growth rate is the greatest for the band, nearest to the height of periodic potential. 
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II. BIPARABOLIC POTENTIAL AND NEAR-TOP BLOCK STATE 

APPROXIMATION 



Stationary Schrodinger equation can be written in a standard form 

[£-,+E-V{z)^^iz) = 0, (1) 

where the particle coordinate z is normalized in units of the potential period L with coeffi- 
cient 27T {z = 27Tz/L = kz, k is the one-dimensional reciprocal lattice constant) and the full 
energy E and the potential energy V{z) are normalized in units of "lattice" recoil energy 
quantum Ej. = h'^k'^/2M, where M is mass of the particle. Taking for convenience the zero 
of the energy on the level of potential energy minimums, biparabolic potential V{z) takes 
the form (see for more details or 0) 

1 -I- (—1 1™ 

Viz) = \ ^ V - (-l)'"x(^ - mn)\ (2) 

where V is the height of the periodic potential, x = 2V"/7r^, m = 0, ±1,±2, and z for 
each m lays in range (m — l/2)7r < z < [m + 1/2)tt. 

The form of the potential is given in fig.[l|, where region / (and others similar to it) will 
be named well-type, and region // (and similar ones) - barrier-type. 

Linearly independent solutions of Schrodinger Eq.(|I]) for well- type region / are written 
in form [0], []5[ 

^i(^) = exp (-^) $ («, ^; V^4) , (3) 

Mz) = exp (-^) * (« + I v^^?) , (4) 
where Zi = z — n, 7r/2 < z < 37r/2, $(..,..;..) is the confluent hypergeometric function and 

The total wave function in this region will be 
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*/(^) = Ci(pi{z) + C2(fi2{z), (6) 

where Ci and C2 are unknown constant coefficients. 

The corresponding solutions for barrier-type region II have the form 



iVx4 \ ^ 1 



lp,{z) = exp y-^^ J * (^^' 2' j , (7) 

= exp (^-yp-^ * (/^ + ^4 -iVxzl) ' (8) 
where 2:2 = — 27r, 37r/2 < 2; < 57r/2 and 

and the total wave function in this region will be 

^Il{z) = Ci^i(^) + C2^2(-2), (10) 

with yet arbitrary coefficients ci and C2. 

With the help of continuity requirements at boundary points z — Sir/ 2 and 57r/2 and 
Bloch periodicity we get the dispersion relation, which can be written in form 

cos(27rP) = 1 + 2Gn{E)G22{E), (Ua) 

or 

cos(27rP) = -1 + 2Gi2(^)G2i(^), (lib) 

where P — p/2hk is the normalized momentum of particle and 

G,j{E) ^ U{z)Tp/{z) + <pi{z)<p^{z)] i,j = l,2. (12) 

where the prime over the function means derivative with respect to z. 

The same requirements, as is well known, with additional normalizing condition 

-3^/2 /•57r/2 

\^j{z)fdz+ / 

'tt/2 J3'k/2 



/■37r/2 /•57r/2 

/ \^i{z)\Uz+ / \^n{z)\^dz = l, (13) 
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define the coefficients ci_2 and ci,2- 

The edge values of energy bands according to ( |llaD and ( |llbD are determined as solutions 
of transcendental equations 

Gi,iE)=0. (14) 

Analysis of these conditions shows that left- and right-side edge points on the positive- 
side axis of quasimomentum (expanded conception of energy bands) are determined from 
conditions Gu{E) = and G 12(E) = respectively for n being even, and by conditions 
G2i{E) = and G22{E) = for n being odd. For keeping analogy with the harmonic 
potential the band with the smallest energy is designated n = 0. 

For even-numbered bands it is appropriate to express the coefficients by Ci_ 

GnjE) + 1 

- G2r{E) e^-P - l""'' 

where Ci is determined from the normalizing condition. For odd-numbered bands it is 
appropriate to express the coefficients by C2; 

_ G2iiE) 

- g;;(e)^^^^^'''' ^^^^ 

= 2GAE)''^ = ^G;M"'' ^''^ 

and use the normalizing condition to determine C2. 

Let's now consider the states, with energies near to the height V of periodic potential: 
E ^ V.ln fig.0 these energies lay inside the dotted lines. For them the parameter /3 (see (^), 
which determines the character of Bloch wave functions in barrier-type region, possesses the 
value (3 = 1/4. The corresponding formulas (|^) and (^ are simplified with help of a well 
known representation for Bessel functions Ju{x) by confluent hj^ergeometric ones: 
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where r(i/ + l) is the gamma function. To use this formula for (|^ we must choose v = —1/4, 
and for - u = 1/A. After corresponding substitutions we arrive to 



To avoid misunderstandings we note, that though hnearly independent solutions ^1(2;) and 
^2(^) for the considered approximation do not depend on energy of the particle, the total 
Bloch wave function "^jjIz) is energy dependant via the coefficients Ci and C2- 

As to wave-functions ipi{z) and (p2{z) of well-type region, we will use the second Tricomi 
expansion for them, 

e-^/^$(a, a + l;x) = r(a + 1) (xx)"'^/^ Yl ^"(^' ^-+n(2v^), (22) 

n=0 

where x = -^^-^ — a, Ao{>c, X) = 1, yli(x. A) = 0, and will restrict ourselves with the 
first member of expansion, which corresponds to sufficiently deep potentials {>c >- 1). The 
respective Bessel functions are expressed by elementary trigonometric functions and as a 
consequence 

(pi{z) = cos(a/E'z), ^2{z) = -^sm{VEz). (23) 

V E 



All the coefficients and dispersion relation (( |lla| ) or ( |116| )) also undergo sufficient changes. 
The dispersion relation, for example, takes the form 

7r/4 2m [J_i/4(m)J_3/4(m) - Ji/4(m)J3/4(m)] cos(7rv^) - J_i/4(m)x 
cos(27rP) = . < 

sm(7r/4) 1^ J,^,(u^^smi^T^/E)-Au^Udu)Js/,iu)^sm{nVE) 

(24) 

where u = 7r^y^/8. 



III. RESULTS OF NUMERICAL CALCULATIONS 



To display the sought dependance, that is the behavior of Bloch wave functions due 
to changing of the energy in hmits of a single (defined) band, we proceed from the above 
mentioned near-top approximation. With the help of computer simulations we first de- 
termined the depth of potential V, so that the last allowed inner-potential band would be 
found immediately near the tops of periodic potential (for fig.(^, with V = 1.4494,the en- 
ergy of mentioned band lays in limits -Emin = 0.3947, -Emax = 1.4494, and when V = 18.65, 
in limits -Emin = 13,64, -Emax = 18,65). Then with the help of formulas (|^) and (pUD, 
with linearly independent solutions (0), (||) and (^), coefficients (|5]), ( |16|) or (0), (|18|), 
dispersion relation (|2^) and normalizing condition (|13|) , we computed the sought values of 
|\E'7-(2;)|^ and |\E'//(2;)|^ as a function of coordinate z along one spacial period. Repeating 
these calculations for all the energy values of the band, we obtain a (continuous) sequence 
of curves, the juxtaposition of which along the energy lines gives us the seeking behavior. 

In Figs. 2a and 2b we present the results of such calculations for two different depths of 
potential; one (Fig. 2a) for a potential depth, containing two energy bands, and the other one 
(Fig. 2b) containing four energy bands. So, Fig. 2a illustrates the energy-space distribution 
of Bloch states in case of relatively shallow potential. Fig. 2b illustrates the same kind 
distribution for relatively deep potential. Left half of the face axis (z) corresponds to the 
well-type region of the potential, and the right half - to the barrier-type region. The energy 
axis is, as is seen, directed from the reader (energy grows from face to back). 

As anomalous one is regarded the behavior of wave function in right-half side(s) of the 
figure (s), that is the dynamics of Bloch functions in the barrier-type regions (see Fig.l) as 
a function of energy. Obviously is seen that the moving along the energy axis from the 
back face to the front face, which physically means getting deeper into potential, does not 
decrease but, quite the contrary, increases the modulo square of the wave function. That 
merely means that the probability of the particle to be found deeper under the potential 
barrier (staying in the limits of a band) is greater than the probability to be found at smaller 
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depths! This kind of behavior is in contrary to the regularity for the case of deepening into 
the periodic potential over various energy bands, or for the case of a single potential barrier. 
That is why we call the behavior anomalous. 

To avoid possible misunderstanding of the mentioned regularity it will be noted that 
there is not any strangeness or abnormality in the fact, that the wave function in barrier- 
type (right) space region is essentially greater in average than in well-type (left) space region; 
it is simply the echo of overbarrier reflecting, when velocity of the particle decreases and the 
wave function (probability of being) respectively increases, in " underbarrier" states. 

The next question that is of interest, is how does the rate of anomalous rising depend 
on the location of the band relative to the height of periodic potential. For this kind of 
calculations the near-top energy approximation, of course, is not applicable and we had to 
proceed from the precise formulas. A detailed analysis of the regularities obtained for Bloch 
state and comparison with single-barrier case we will give in another publication. We would 
like only to note that while moving away from near-top energies as inside as outside the 
periodic potential, the rate (scaled in energy units) of anomalous growth decreases. 

Finally we note that we carried out the anomalous behavior also for the Kronig-Penny's 
potential, but for that case the anomaly is much weaker, than for the case of biparabolic 
potential (see fig.(^). It is easy to understand since the single barrier transparency in high 
energy region, and hence the effectivity of matter-wave interference is bigger for biparabolic- 
form potential than for rectangular-form one. 

This work was supported by ISTC Grant A - 215 - 99. 
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FIGURES 



FIG. 1. Biparabolic potential and identification of well-type(I) and barrier- type (II) regions. 

FIG. 2. Modulo square of Bloch wave functions on the plane energy (E)-coordinate (z) for a 
near-top (with -Emax = V) energy band. The depth (heightjof potential is F = 1.4494(a) and 
V = 18.65(6), scaled in units of recoil energy Er = (2hkf/2M. The mentioned feature is in detail 
described in the text. 

FIG. 3. Modulo square of Bloch wave functions on the plane energy (E)-coordinate (z) for a 
Kronig-Penny potential. The depth of potential \sV = 1.668, scaled as in Fig.2. 
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